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Abstract. The first step toward the application of an effective non partial wave (PW) numerical approach to 
few-body atomic bound states has been taken. The two-body transition amplitude which appears in the kernel of 
three-dimensional Faddeev-Yakubovsky integral equations is calculated as function of two-body Jacobi momen- 
tum vectors, i.e. as a function of the magnitude of initial and final momentum vectors and the angle between them. 
For numerical calculation the realistic interatomic interactions HFDHE2, HFD-B, LM2M2 and TTY are used. 
The angular and momentum dependence of the fully off-shell transition amplitude is studied at negative energies. 
It has been numerically shown that, similar to the nuclear case, the transition amplitude exhibits a characteristic 
angular behavior in the vicinity of "'He dimer pole. 



Introduction 

In recent years the ^He trimer and tetramer have been the 
center of several theoretical investigations (see, for exam- 
ple, Refs. II1I2I3I and references therein). From all em- 
ployed methods in these studies, the Faddeev-Yakubovsky 
(FY) schemas are perhaps most attractive since they re- 
duces the Schrodinger equation for three (four) particle 
systems into a coupled set (two coupled sets) of integral 
or differential equations which can be used to study the 
bound and scattering states in a rigorous way. The differen- 
tial form of FY equations has been successfully applied in 
nuclear bound states calculations, but there are limitations 
in its application to atomic systems. The limitation arises 
from eccentricities of the interatomic interactions, since in- 
teratomic interactions often contain very strong short range 
repulsion which leads to tedious and cumbersome numer- 
ical procedure. In calculations of atomic systems, because 
of the short range correlations, one needs a large num- 
ber of PWs to obtain the converged results. To overcome 
this problem few numerical techniques are developed, the 
tensor- trick method [TTl, representation of Faddeev equa- 
tions in Cartesian coordinate 1 6 1 , the operator form of Fad- 
deev equations in total angular momentum representation 
Q and also a hybrid method 1 8 1 . The three and four-body 
atomic bound states have been also studied with short-range 
forces and large scattering length at leading order in an Ef- 
fective field theory approach [9|-111J, but these investiga- 
tions are also based on PW decomposition and the interac- 
tions are restricted to only s-wave sector. 

By these considerations we are going to extend a nu- 
merical method, which has been successfully applied to 
nuclear bound and scattering systems and avoids the PW 
representation and its complexity, to atomic bound states. 
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It should be clear that the building blocks to the few-body 
calculations without angular momentum decomposition are 
two-body off-shell transition amplitudes, which depend on 
the magnitudes of the initial and final Jacobi momenta and 
the angle between them. Elster et al. have calculated the 
A^A^ transition amplitude for spinless particles in a non PW 
representation by using the Malfliet-Tjon type potentials 
lfT2l . Our aim in this paper is to calculate the matrix ele- 
ments of the fully off-shell two-body transition amplitude 
at negative energies for realistic interatomic interactions, 
we study the momentum and angle dependence of transi- 
tion amplitudes. This paper is organized as follows. In sec- 
tion[T]we represent the explicit form of studied interatomic 
interactions in configuration and momentum spaces. In sec- 
tions|2]and[3]we present our numerical results for homoge- 
nous and inhomogenous Lippmann-Schwinger equations 
in a non PW representation. An outlook is provided in sec- 
tion IH 



1 '^He-'^He Interatomic Interactions 

In this study as ^He-^He interatomic interactions we use 
the realistic HFDHE2 fill, HFD-B [T4l, LM2M2 US and 
TTY 1161 potentials. The semi-empiiical HFDHE2, HFD- 
B and LM2M2 potentials, which are constructed by Aziz 
and collaborators, have the general form 



V(r)^s\Va{x) + Vb(x)l 



(1) 



where x - —, r and r,„ are expressed in the unit A. The 
terms Va{x) and Vbix) read 



'An sin 



Vaix) = 



Inix — xi) n 
X2 — xi 2 

0, 



-Hi , Xi < X < X2 



(2) 



X i {Xl,X2\ . 
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Fig. 1. The interatomic "'He-'^He potentials as a function of ttie distance r between the atoms. The region around minima of the potentials 
is shown in the inset of the figure. 



C(, Cg c 



and the function F(x) is given by 



F(x), 



Fix) 



(3) 



(4) 



1, 



x> D. 



The parameters of the HFDHE2, HFD-B and LM2M2 po- 
tentials are given in Table [T] The explicit form of the theo- 
retical TTY potential is 



V(r) = A (yex(r) + VdispCr)), 



(5) 



where r stands for the distance between the "^He atoms 
given in atomic length units. The function VexCO has the 
form 



7 

with p - 1 • The function Vdispir) reads 

zp 

N 

Vdisp(r) = -2C2„/2„(r)r-2«, 

«=3 



(6) 



(7) 



the coefficients C2,, are calculated via the recurrency rela- 
tion 



L2«- 



(8) 



and the functions /2„(r) are given by 
where 

b(r) ^2/3- 



k=0 



(9) 



(10) 



The parameters of the TTY potential are given in Table |2] 
In Figure ([TJ all potentials are plotted as function of dis- 
tance between the "^He atoms. Moreover, in the inset to this 
figure the region around the minima of the potentials is 
shown for providing a better comparison. 

In order to be able to implement the introduced inter- 
atomic interactions in few-body atomic bound and scat- 
tering state calculations in momentum space, we need to 
transform these potentials to momentum space. The matrix 
elements of the potentials can be obtained by following re- 
lation 

Vip,p') = V(p,p',Xpp,) 
I f"' 

^r- drrsm(qr)V(r) ; ^ = |p - p | 
^^q Jo 



2^ Jo 



drr^V(r) 



(11) 



where p and p' are magnitudes of initial and final two-body 
Jacobi momentum vectors, Xpp' is the angle between them 
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Fig. 2. Momentum-space representation of interatomic ''He-''He potentials V(p, p',Xpp>) in fixed angles Xp^t = 0, ±1. 



and q - (p^ + p'^ - Ipp'xpp')^ is the difference between that the potentials vanish at enough large values of Jacobi 



them. Clearly the knowledge of the structure and range of 
potentials is important in the calculation of two-body tran- 
sition amplitudes, which appear in the kernel of the integral 
equations of two-, three- and four-body bound and scatter- 
ing calculations. In Figure (|2]i the momentum dependence 
of the potentials are shown at fixed angles Xpp' - 0, +1. 
As shown all potentials have similar behavior The ridge 
around p = p' arises from strong repulsive core. The be- 
havior of the potentials at forward angle, i.e. x 



pp 



+ 1, is 

different from other angles, and according to Eq. ( fTTb the 
value of the potentials in this angle is fixed for p - p' . Note 



momenta, i.e. p 



= 20 A . We should mention that for 



calculation of matrix elements of potential, Eq. (fTTI) . the 
cutoff value of r'™^ = 30 A and 200 mesh points have been 
used. 
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Table 1. The parameters of the ''He-^'He interactions which are 
constructed by Aziz and collaborators. 



Parameter 


HFDHE2 


HFD-B 


LM2M2 


£[K] 


10.8 


10.948 


10.97 


r„, [A] 


2.9673 


2.963 


2.9695 


A 


544850.4 


184431.01 


189635.353 


a 


13.353384 


10.43329537 


10.70203539 


P 





-2.27965105 


-1.90740649 


Ce 


1.3732412 


1.36745214 


1.34687065 


Cs 


0.4253785 


0.42123807 


0.41308398 


Cio 


0.178100 


0.17473318 


0.17060159 


D 


1.241314 


1.4826 


1.4088 


Aa 






0.0026 


Xi 






1.003535949 


X2 






1.454790369 



Table 2. The parameters of the ^He-^He TTY potential. 



A[K] 


315766.2067 


/3 [(a.u.r'] 


1.3443 


D 


7.449 


N 


12 


Cfi 


1.461 


Cg 


14.11 


Cio 


183.5 



2 "^He Dimer 

The ^He dimer can be described by homogeneous Lippmann- 
Schwinger integral equation: 



-J 



(12) 



This integral equation can be solved numerically by direct 
or iterative methods. We have solved this integral equation 
by direct method and the numerical results for dimer bind- 
ing energy by using the introduced interatomic interactions 
are given in Table (O in comparison to corresponding PW 
and experimental results. In our numerical calculations the 
"^He atom mass is defined by — = 12.12KA^. For dis- 
cretization of the continuous momentum and angle vari- 
ables we have used the quadrature Gauss-Legendre by us- 
ing linear mapping for all variables. The number of mesh 
grids for Jacobi momenta and angle variable are 200 and 
150 correspondingly. 



Table 3. The calculated dimer binding energy in unit of mK for 
realistic interatomic potentials in comparison to corresponding 
PW and experimental data. 





HFDHE2 


HFD-B 


LM2M2 


TTY 


PW|17| 
Presnet 


-0.83012 
-0.83011 


-1.68541 
-1.68540 


-1.30348 
-1.30347 


-1.30962 
-1.30962 


EXP |T8| 




-1 1+*" 

'■'-0.2 







3 Two-Body Transition Amplitude 

The building blocks for few-body bound and scattering 
state calculations are two-body transition amplitudes T which 
follow the inhomogenous Lippmann-Schwinger equation 



T + VGoT, 



(13) 



where V is the two-body, e.g. two-atom, potential and Go = 
(E - //()) ' is free two-body propagator. For momentum 
space calculations one needs the matrix elements of transi- 
tion amplitude in desired energy E which can be obtained 
by representation of Eq. (fT3]) in two-body basis states l,12J 



r(p',p;£) = y(p 



',p) + J d'p" 



-;rT(9 ,p;£)- 



(14) 



In order to solve this three-dimensional integral equa- 
tion directly without employing PW projection, we have to 
define a suitable coordinate system. To this aim we choose 
vector p parallel to z-axis and vector p' in the x - z plane 
and express the integration vector p" with respect to them. 
By this considerations Eq. ( fT4l l can be written explicitly as 



Tip' ,p,Xppr,E) = V{p',p,Xpp,) 

+ dp"p"^ dxpp„ dip" —V{p\p",x, 

Jo J-i Jo £ - £- 

m 

xT{p",p,Xpp„;E), 
where 

P P". 



p' p") 
(15) 



Xp' p" 



P ■ P — Xpp'Xpp' 



(16) 



The (fi" integration acts only on V{p', p" , Xp' p"), so this 
integration can be carried out separately as 

V{p' ,p" ,Xpp>,Xppn) = I dip"V{p' ,p" ,Xp> pn), (17) 

Jo 

and consequently the integral equation (fTSl l can be written 



T(p',p,Xpp>;E) ^ —v{p\p,Xpp.,\) 
Zn 



dxp p" 



Jr»oo pi 

xT(p",p,Xpp,r,E) 



v(p',p",Xpp>,Xpp") 



E-^ 



(18) 



For a specific value of the off-shell momentum p and 
energy E and after discretization of continuous momentum 
and angle variables, this two-dimensional integral equation 
can be turned into a system of linear equations as AT = B, 
where A and B are composed of kernel of integral equation 
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and potential matrix elements respectively. For our numer- 
ical calculations we use the Lapack Fortran library lfT9l 
to solve the obtained system of linear equations. Certainly 
for few-atomic scattering state studies one needs to cal- 
culate the transition amplitude at positive energies which 
leads to a singularity in free propagator This Cauchy sin- 
gularity can be splitted easily into a principal-value inte- 
gral and a 5-function imaginary part. Since we are going to 
use the recently developed formalism for three- and four- 
body bound states f20l-f23l for "^He trimer and tetramer 
calculations we study the behavior of transition amplitude 
at negative energies. To this aim we have solved the two- 
dimensional integral equation (fTSl l by using 40 and 41 mesh 
grids for Jacobi momenta and spherical angles variables re- 
spectively. We would like to mention that by considering 
the symmetry property of the angle argument Xp' p", the 
polar angle integration in Eq. ( fTTl i can be done on inter- 
val [0, |] by using 10 mesh grids. Our numerical results 
for fully ofFshell transition amplitudes T{p, p' , Xppt\E) are 
shown in Figure Q at energy E - -IQQniK, which is 
close to "'He trimer binding energy, in fixed angles Xppi - 
0, ±1. In Figure (|4]i we have shown the momentum and 
angular dependence of half offshell transition amplitude 
T{p, po - ylm\E\, x\ E) at energy E = - 100 inK. 

As shown in Ref. lfT2l the bound states of two-body 
system lead to poles in the transition amplitude and the an- 
gular dependence of transition amplitude exhibits a very 
characteristic behavior in the vicinity of the bound state 
poles, which is given by the Legendre function correspond- 
ing to angular quantum number of the bound state 

r(p,p;£) — Piixpp,)— — - — , (19) 

47r E - Eh 

where 

g,(p)^ dp'p'^v,(p,p')Hp'), (20) 
Jo 

where vi(p, p') and i//iip') are partial wave components of 
potential and tw-body wave function. To investigate this 
characteristic behavior in atomic case we have shown in 
Figure (|5]l the angular dependence of transition amplitude 
Tipo, Po, x; E) at the energy range -400 <£■<-! niK. 
Clearly we can see by decreasing the magnitude of energy 
the angular behavior of transition amplitude is correspond- 
ing to the zeroth Legendre polynomial, i.e. Po{xpp>). For 
a better representation of this angular behavior we have 
shown in Figure (|6]l the transition amplitude just for few 
energies close to dimer i-wave pole and also the magni- 
tudes of transition amplitude at energy £ = - 1 mK are 
Usted in Table Q. 



4 Outlook 

The first step toward the application of an established non 
partial wave approach to few-body atomic bound states has 
been taken. The necessity of using this non partial wave ap- 
proach comes from this fact that in few-body atomic cal- 
culations one needs a large number of partial wave com- 



Table 4. The matrix elements of T{po, po, Xp,/ ; E) in unit of A'.A' 
with pq = y/m\E\ at energy E = - I mK, which is close to .s-wave 
"^He dimer pole. 

x„„, = p.p' HFDHE2 HFD-B LM2M2 TTY 



ponents, which is caused by very strong short range re- 
pulsion of interatomic interactions, to reach proper con- 
verged results. Instead of using standard partial wave rep- 
resentation which leads to tedious and cumbersome numer- 
ical procedure we intend to extend a non partial wave ap- 
proach which has been successfully applied to few-body 
nuclear systems. In the first step toward this goal the ma- 
trix elements of transition amplitudes which appear explic- 
itly in the few-body calculations have been calculated di- 
rectly as function of two-body Jacobi momentum vectors. 
The calculated matrix elements can be entered in kernel of 
three-dimensional Faddeev- Yakubovsky integral equations 
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-0.62977 


-2.5366 


-3.0905 


-3.0486 


-4.8682 
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-0.43872 


-2.5373 


-3.0912 


-3.0493 


-4.8689 
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Fig. 3. Momentum dependence of the fully ofFshell transition amplitude T(p, p', Xppi ; £) at £ = - 100 mK in fixed angles Xppi = 0, ± 1 . 



to study the "^He trimer and tetramer ground and exited 
states. The numerical calculations for these bound states 
are currently underway. 
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